Skip to main content

Full text of "Convergence of the self-dual Ginzburg-Landau gradient flow"

See other formats

Convergence of the self-dual Ginzburg- Landau gradient flow 

Sophia Demoulini 

Centre for Mathematical Sciences, Wilberforce Road, Cambridge, CBS OWB, England 

email:sd290@cam. ac. uk 


We prove convergence of the gradient flow of the Ginzburg-Landau energy functional on 
a Riemann surface in the self-dual Bogomolny case, in Coulomb gauge. The proof is direct 
and makes use of the associated nonlinear first order differential operators (the Bogomolny 
operators). One aim is to illustrate that the Bogomolny structure, which is known to be of 
great utility in the static elliptic case, can also be used effectively in evolution problems. We 
also identify the minimizers and minimum value of the energy when the Bogomolny bound is 
not achieved (below the Bradlow limit). 

MSG classification: 58J35, 35Q56 

1 Introduction and statement of results 

We consider the self-dual Ginzburg-Landau, or abelian Higgs, energy defined over a two dimensional 
Riemann surface S: 

V(A, = - / I^A^ A *Dx^ + F*F+-{l- |$|^))^ * 1 

Here the Higgs field, is a section of a complex line bundle L — )■ S with fixed metric h, and 
F>A = (-Ci, D2) = V — iA is the covariant derivative operator defining an connection on L, with 
curvature 2-form —iF. We will assume that the degree N = deg L is a non-negative integer (but 
everything in this article has an analogous version in the negative case). It is known from [2] that 
for |S| > 4ttN the minimum value of V is ttN and is attained on the set of self-dual vortices, which 
are solutions of the Bogomolny equations (11. 5p ; these form a system of first order partial differential 
equations, solutions of which are minimizers of V when they exist. The presence of such a system of 
first order equations is a special feature of the functional V which is related to self-duality ([3 [2]). 
The main aim of this article is to show how it is possible to prove convergence of the corresponding 
gradient fiow quite simply by making use of the Bogomolny (or self-dual) structure in (jl.4p - (jl.5p 

it is also proved that for | S | < AttN the minimum value of V is ttN + where I = 27riV — '-^ , 
and this is achieved by taking <I> = and A to be a constant curvature connection. 

A similar proof of convergence for the SU{2) Yang-Mills-Higgs functional on was carried 
out in [7]; this is closely related to the case |S| > 4ttN in this paper. The crucial structural 
features exploited there for Yang-Mills-Higgs also hold for the Ginzburg-Landau flow in the case 
|S| > AttN: by (j3.18p the first order Bogomolny operators (j3.16p - (j3.17p converge to zero, while 
()3.30p then implies the convergence of (A, $) as t — t- -|-c«. The case |S| < iirN is different: indeed 
this latter inequality is an obstruction to the existence of solutions to the Bogomolny equations 

([2]), and so the first stage of the argument just outhned necessarily fails. We give, however, an 
alternative set of quantities in (I3.45p . whose large time behaviour can be analyzed, and which serve 
as an effective alternative to the Bogomolny operators in this case. Compared to the adiabatic 
approximation method in [4j, where convergence of the Ginzburg-Landau flow was proved for the 
case S = M^, the proof here is certainly more efficient although it makes less contact with the 
physics of vortices. It is also possible to prove convergence by the very general method based on 
the Lojasiewicz inequality in as has been done in [6j. In comparison with these previous results, 
the present method does provide the additional information that convergence is at an exponential 
rate. Furthermore in view of the importance of the Bogomolny self-dual structure in the elliptic 
case, it seems worthwhile to illustrate its utility in the parabolic context, as was done also in [3] 
for vortex dynamics in a conservative (Hamiltonian) context. 

Introducing conformal coordinates {x^}^^^, in which the metric takes the form 

g = gjkdxHx^ = e^P{{dx^f + {dx^f) 
with associated area form d^g = e^^dx^ A dx'^, the energy functional V is given by 

V(A,cI>) = i|^(|DA^|V2p + B'- + \{l-\<^\^f) dt^g. (1.1) 

If we fix a smooth connection V on L then Da is determined uniquely by a real 1-forni 

A = Aidx^ + A2dx'^ G n\^) 

where : S — )• R, according to 

I?A:=V-iA. (1.2) 

The curvature, or magnetic field B, of a connection Da = Djdx^ is determined by DjDj^^dx^ A 
dx^ = —iF^ = —iB^dfig. Its integral is a topological invariant of L: 

I Bdjjg = 2t:N (1.3) 

with = deg L. We will choose V to have constant curvature 6; this value is then fixed topologically 
as 6 = 27rA^/|S|, where |S| = dfig is the area of S. It follows that 

B = b + *dA = b + e-^P{diA2 - ^2^1). 
1.1 The Bogomolny equations 

A crucial property of the energy functional (jl.ip is that (in appropriate functions spaces in which 
e.g. integration by parts is valid) it admits a Bogomolny decomposition into a sum of squares of 
first order terms: 

V(A,$) = ^£(4|9A$pe-2^ + (S-l(l-|$|2))2) dt,g + ttN (1.4) 

where as above N = degL and 

dA<^> = ^iDi+iD2)'^. 

If the following first order equations, called the Bogomolny equations, 

Sa* = 0, 

have solutions in a given class, they will automatically minimize V within that class (this theory - 
and conditions under which solutions exist - is developed for the plane and for surfaces in [8lll01[2]). 

Our main theorem concerns the large time behaviour of global weak (Sobolev) solutions to 
1.61 in Coulomb gauge (the existence and uniqueness of which is shown in section [2]). We show 
that as t — 7- oo this solution converges exponentially fast to a minimizer of (jl.ip a solution of the 
Bogomolny equations when |S| > 4ttN . We make the assumption of small (close to equilibrium) 
initial energy. 

We will make use of the Bogomolny equations to derive the asymptotic convergence of weak 
solutions to the gradient flow equations in section [3l In fact, the convergence of such solutions to 
minimizers of the energy provides an independent proof of the existence of minimizers and solutions 
to the Bogomolny equations in the case when |S| > AttN. 

1.2 The time-dependent equations and statement of the main theorem 

Associated with the energy functional (jl.ip are the Ginzburg-Landau gradient flow equations in 
the variables (A, = (^o,A,^>) = (^0,^1,^2,^) on [0, 00) x S: 

(|,A)(A,cI>) = -V(A,cI>) 

and equivalently, 










- Aa$- 


where the covariant Laplacian Aa is deflned as 

This is the fully gauge invariant version of the Ginzburg-Landau gradient flow. The Higgs field ^, 
now depending on time t > , is a section of a complex line bundle L = [0, 00) x L — t- [0, 00) x S 
on which there is a connection 

A = (^0,^1,^2) = (^o,A) 

with associated covariant derivative D = {Dq, Di, D2) = {Dq,Da), where Dq = §1 — iA^ is the 
time component and Dx = {Di,D2) is the spatial component as in (jl.2p with Ai now depending 
on time t > as well as x G S. The dependent variables are thus the Higgs field $ and a real 
1-form AQdt-\- Ajdx^ on [0, c«) x S. The system is invariant under the infinite dimensional group of 
gauge transformations: if xi^i^) is a smooth function then [Aq, Ai, A2,^) is a smooth solution of 
(jl.6p if and only if {Aq + dtX-, Ai + dix, A2 + d2Xi e*^<I>) is a smooth solution of (|1.6p . This means it 

is possible to place further conditions on the solution, and this is necessary to obtain uniqueness. 
We impose the Coulomb gauge condition : 

div A = e-^P{diAi + ^2^2) = ; (1.7) 

a solution, or set of initial data, is said to be in Coulomb gauge when this condition holds for all 
relevant t. In addition, we will require that 

Aodfig = (1.8) 

at all relevant t: this may be achieved by applying to the solution {Aq, Ai, A2, $) at each (t, x) the 
gauge transformation x{t) = fT,^odfJ.g- 

Initial data and Coulomb gauge: We specify (A, <1>) = {Ai,A2,^) at t = which in Coulomb 
gauge also determines Aq initially. To see this, we eliminate Aq from the full (gauge invariant) 
equations (11. 6p by taking the divergence in the equation for A under the condition (11. 7p . This 
leads to the equation for Aq in Coulomb gauge, 

- AAq = div Da<^) . (1.9) 

Lemma 1.1 For each fixed t > there is a smooth function (A, <I>) 1— )• a(A, <I>) from H^{L(B^^) — )■ 
H^'*^, for any 1 < q < 2, which is the unique solution of (\1.9h satisfying (II. Sp . 

Proof Because ^ e then for all p < 00, <I> G LP, Da^ G L"^ and so div Da^) G H~^''^ for 
all g < 2. It follows that (|1.9p has a unique solution = a(A, <I>) verifying (jl.Sp as described in 
the lemma. □ 

Remark 1.2 Furthermore, if ^ satisfies (jl.6p then fort a.e. —AAq = div {i^, Da^) = {i^,DQ^) , 
and so by standard elliptic theory, 

\\Ao\\h2,p < Cs\\{i^,DQ^)\\LP , 

where finiteness of the right hand side will be justified in section Wl\ By Sobolev imbeddings, here 
applied for 1 < p < 2 at each time t > 0, 

WAqWh^.p < CiW^W 2p ||Do^||l2 and so \\dAQ\\i2 < C2\\^\\Hi Wo^Wl^- (1-10) 

L-'-P ^ 

As a particular consequence, ^0 is initially determined through the initial data for (A,<I>). In 
addition, the estimates p.lUp will be used below to derive estimates for Dq^. 

Explicitly applying the gauge conditions (jl.Tp and (jl.Sp and substituting Aq = a(A, $), the 
system (jl.6p becomes the following nonlocal parabolic system: 

A = AA + da + {i^, Da^) 

$ = A$ - 2ie~^Pj2^j^j'^ - e"^''| A|2$ + + ^(1 - |$|2)^. (l'^^) 

with initial data specified for (A,$). We discuss existence of this system in section [2] and we will 
need the following spaces. 

Function spaces: We will work with spaces of functions H'^'P(T,) which may be defined either 

(i) by using a partition of unity to reduce to the Euclidean case (the standard spaces H^'P{M?)), or 

(ii) by appealing to the spectral theory for —A to define a functional calculus and then introducing 
the norm as 

||/||H^,P = ||(l-A)f/||iP. 

Let fj be an orthonormal basis of eigenfunctions, —^fj = ^jfj with < Ai < A2 < • • • , then for 
p = 2 an equivalent norm on H'^ = H'^'^ is 

11/111,.= ^(1 + A,r|c,f, where / = ^c,/,. (1.12) 


These definitions extend to sections of vector bundles L, in which case we write H'^{L) for the 
corresponding spaces defined with respect to the connection V and its corresponding Laplacian 
g-2p^Y2 _|_ y2-j. £qj. 5 = X an equivalent norm is 


\\nHiiL) = jj\'^\' + \^mdfig. (1.13) 

In the above integral the inner products are the standard ones induced from h and g. The definitions 
also extend to 1-forms using the Hodge Laplacian A = —{d5 + 5d), (with the sign chosen so that 
it is non-positive). In conformal coordinates: 

{AA)i = diie-^PidiAi + d2A2)) - eijd,ie-^PidiA2 - Qa^i)) 
= -eijdjB 

when the Coulomb gauge (jl.7p holds. The corresponding norms are written H''^{Q^) where is 
the space of one- forms on S with coefficients in . 

The main theorem 

For data we derive existence of a solution also in in the following section in the Coulomb 
gauge. The main result concerns the time asymptotic behaviour of these solutions which are shown 
to be exponentially converging to an energy minimizer. The minimum energy depends on whether 
the size of E is large enough to support the existence of vortices (O [TO]): 

■nN if S > A-kN 

^^^(4.iV |E|)^ .f^<4.Ar. 



Vo = V(A(0),$(0)) 
be the initial energy. The following is the main theorem 

Theorem 1.3 (Main Theorem) Given finite energy initial data (A(0),$(0)) G H^{Q}®L) sat- 
isfying the Coulomb gauge condition (jl.7p and such that ep with 

< Vo — Vmin = eo 

is sufficiently small there exists a unique global solution (^0, A,^>) with Aq € C([0, 00)); i?i'P(S)) 
for any p < 2, and (A, ^>) G C([0, 00); ^^(f]^ L)) of ([II]), satisfying ([LTD and ([TS]) with the 
following large time behaviour: 

> iirN then {A, ^ {Aoo,^oo) in {n^e L) where (Aoc^oo) is a solution of the 
minimum energy static (Bogomolny) equations and where <I>oo has precisely N zeros (vortices) 
on E. Also ^o(^) ~^ strongly in L"? for all q < oo. The minimum value of the energy in 
H^{L is attained and is ttN. 

i/|E| < A-kN then{A,^) ^ {A^,Q) whereD = V-iA oc is a constant curvature connection 
at which the minimum value of the energy in H^{L © O-*^) is achieved and equals 

V(A„.0)= min V-^A' + a^i^. 

m{L(Bn^) 8|S| 

Convergence is at an exponential rate, namely, 3 c,6>0 depending only the initial data such that 

||(A,$)-(Aoo,$oo)||l2 <ce~^*, respectively, || (A, $)- (Aoo, 0)11^2 < c e-"^*. (1.14) 

Global existence for (jl.6p is established in the following section, making use of standard (semi- 
group) techniques (theorem 12.10 . The proof of asymptotic convergence, which is the main content 
of theorem 11.31 is then given in section [3l 

Remark 1.4 Convergence in stronger norms for s < 1 can be deduced by interpolation between 
LP' and from (jl.l4p . using convergence in ("low" norm) and boundedness in (or "high" 
norm) guaranteed the energy non-increase. Furthermore, when boundedness can be proved in higher 
norms than then interpolation would lead to convergence in for higher s by standard parabolic 

Remark 1.5 The case 

|S| = 47riV 

is degenerate (as becomes clear from reading the proof) and it is conceivable that exponential con- 
vergence may not occur, even if the initial energy is close to the minimum energy. 

Remark 1.6 The condition that the initial energy be close to the minimum energy can be relaxed 
in situations when there are no non-minimal critical points (JBi chapter 3]). The energy is strictly 
decreasing (as it is a gradient flow) and so it must decrease to the energy of a critical point. If non- 
minimal critical points do not exist then the energy must decrease to its minimum value, in which 
case it approaches the minimum energy, and then the conclusion of the theorem holds regarding 
exponential convergence. This is not discussed here further. 

2 Existence theorem 

There are different ways to obtain a global existence theorem for these equations, for example using 
maximum principles (as for the case of the same equations on in [4]) or energy methods as done 
here, using semi-group techniques explained in e.g. [12] and elsewhere. To use the energy norm 
method, we first reduce our system in Coulomb gauge to a system for (A, <I>) given in (11.110 : then 
we show that for data (A(0),$(0)) there is a unique global weak solution of those equations, 
also in H^, varying continuously in time. 

Theorem 2.1 (Existence in Coulomb gauge) Given initial data (A(0), <I>(0)) G H^{Q^ © L) 
(which implies finite initial energy V(A(0), <I>(0)) = Vq < oo), there exists a unique mild solution 
of (jl.lip with regularity (A,<I>) G C([0, oo); © L), which satisfies 


for all t > 0. For smooth initial data the solution is smooth, and for H^{U^ © L) initial data the 
solution is the limit in C;oc([0, oo); H^{Q,^ © L) of smooth solutions. 

Proof Writing U = (A, $) and CU = (A A, A$) the system (fTTTI) is of the form U = CU + FiU), 
which can be treated as an abstract evolution equation in the space X = H^{Q}'(BL). We introduce 
also the auxiliary space Y = H~^[Q^ © L), and make use of the following facts: 

• e*"^ is a strongly continuous semi- group of contractions on X for t > 0, 

• e*^ : y ^ X for t > 0, with ||e*-^||y_,x < Cr^/^ for < t < 1, 

• T : X ^ Y is a smooth function satisfying \\T{Ui) — J^{U2)\\y < K{R)\\Ui — U2\\x for 

The first of these is a standard property of the heat equation. The second can be derived using the 
spectral representation introduced above: 

which implies that 

\\e'''f\\ls=Y,e-''^^{l + \,r\c,\''<Cs,rmffHr 


with Cs,r{t) = sup^>Q — ~ ~ as t — )• 0. The third assertion follows by examining 

the various terms which constitute and applying appropriate embeddings for the spaces, in 
particular C and Ls c H^^ (continuous embeddings). For example, the current term 
(i<I>, V<I>) in the $ equation arises as a continuous bilinear map x L'^ ^ La c H^^ bounded by 
|(i<I>,V<I>)| 4 < c||<I>||^4 ||V<5||j;^2. The other terms in are treated similarly (using lemma [TTT] to 
handle the a terms). 

These three properties imply (see [121 §15.1] that the integral operator U — e^^^^''^J^(U{s))ds 
is a contraction on C([0,T];X) for sufficiently small T, and in fact T can be taken to be a 
positive non-increasing function of ||Z^(0)||x- This implies that there is a unique local solution 
U G C{[0,T]; X), and also that, given two sets of initial data Z^i (0) , (0) in X, the well-posedness 

mayijUi{t)-U2{t)\\x < c\\Ui{0) - U2m\x 

holds for some c > 0, for T sufficiently small (again depending only on ||ZYj(0)||x-) From this it 
can be deduced that for initial data the solution remains in H^, and in fact for initial 
data the solution remains in H'^ for all s > 2. Thus the solutions can be approximated by 
regular solutions, and obey the energy non-increase : V((A, <I>)(t)) < V((A, <l>)(s)) for t > s. 

As a consequence the norm is globally bounded and there exists a unique global solution 
U G C([0,cx));X) of the corresponding integral equation 

U{t) = e'^U{0) + f e^'-'^'^F{U{s))ds 

with the approximation property asserted. □ 

3 Asymptotic behaviour and proof of the main theorem 

We continue to use the Coulomb gauge in which Aq is eliminated as an independent variable by 
lemma [Tm We consider the cases |S| > iirN and |S| < AttN separately. In each case we introduce 
auxiliary variables which are specially tailored to reveal the asymptotic convergence of (A, to 
a limit {Aao,^oo) which is characterized differently depending on whether the surface area |S| is 
bigger or smaller than AttN. In the following proofs we shall make use of the norm 

= [ (|7?|2 + \Dj^rj\^) dfig (3.15) 
defined with respect to the time dependent connection = V — iA. 
3.1 Proof in the case |S| > inN 

It is useful to introduce the Bogomolny variables v = v{A, <I>) and ij = r]{A, defined as 

r] = Ba^ (3.16) 
v = B-^il-\^\') (3.17) 

where, as above, Ba = ^{Di + iD2). In terms of these variables the energy decomposes: 

V = ^JjA\rj\h-^P + v^)df,g + ttN, 

as was seen above in ()1.4p . The evolution of these variables is according to 

DqT] - 4:8 Aie'^^d All) + = -vr], 

d (3-18) 
(|-A + |cI>|> = -4H^ 

These equations can be obtained from (jl.6p by first applying the operator Ba to the equation for 
^, using the identities (in conformal co-ordinates) 

48aBa = Dl + Dl + e^PB , Aa = e-^f'iDl + DI) 


9aAa - 4BA{e-''PdA)BA = -{BB) - BBa 



(where BAq = {BAq) + AqB) and finally using the equation for A in the last commutator. (To 
derive the final commutator use also that the background connection V is fixed independent of t, 
i.e., [|,V] = 0.) 

We will show (using lemma [3^ that 

(??, v) — > (0, 0) as t — )• oo 

(in specified spaces) and then deduce via lemma [3l4l the convergence of (A, as stated in the main 
theorem. For this consider two energy-type quadratic forms, respectively associated to each of the 
equations in 13.181 (at each time t > 0) defined as 

QaAv) = j^{4.e-^'\dAV\^ + e-^''\'^\^\v\^)di^g , (3.19) 
Q'fiv) = j (^e'^P{\Div\'^ + \D2v\'^) + \<l>\'^v^yfig. (3.20) 

Remark 3.1 It is shown in ^3] [lemma 3.2.2] that these quadratic forms are coercive in and 
respectively, provided that > m > 0.' to be precise, under this assumption on $, there 

exists 7 = 7(m, M) > (where recall that V < M < oo and M = ttN + 1), such that 

Qa, $(??)> 7||??||^^ and Q<i>('u) > 7||v||^i • (3.21) 

That the assumption on \\^\\l2 is a valid one can be shown as follows: by the Bogomolny decompo- 
sition of the functional V (II. 4p and the topological invariant (11. 3p . 

ISI 11$ 

27riV - — + 



from which it follows that. 

2 2 

\q>\\\^ = |S| -AttN + 2 I V 

> |S| -47riV-2(2eo)^|S|^ ^^'^^^ 


for eo sufficiently small. 

In terms of the quadratic forms Qa,^, the variables rj, v satisfy the following inequalities: 

Lemma 3.2 Let (A, <I>) be as in theorem \2.1\ and define {rj,v) as in p.l6p - (j3.17p . Then for any 
time T > 0, the following identities hold 

e''^UT)\\l,+2 f e^'\QAAv)-m\h)dt < ||77(0)||i. - 2 j"" J^e'''v\v\'e-'Pdi,gdt, 



e''^''\HT)\\l2+2 I e^^'{Q^{v)-6\\v\\l2)dt < \\v{0)\\l2 - 8 j l^e^^'v\ri\^e-^PdfXgdt. 


Proof For smooth solutions multiply 13.181 respectively by t] and v and integrate with respect to 
e'^^^e~^Pdfigdt over [0, T] x S. For more general finite energy solutions use the approximation prop- 
erty in theorem 12.11 □ 

The above lemma implies the exponential decay of r], v under the conditions on $ for co-ercivity 
and for initial energy sufficiently close to its minimum. 

Corollary 3.3 (Exponential decay) For all 6 e (0, ^) and < eo with ^/2eo < g^, ifV{A, <!>)- 
ttN < eo, then 


snpe''\Ut)\\l + \\v{t)\\l) + 6 e''\Ml.+\\vfm)dt< 2(||r?(0)||i. + ||^(0)||i.) . (3.25) 
f>o Jo ^ 

Proof Recombining the terms in ()3.23p we have the following bounds 

[QaMv) - SMh) dt + 2^ 1^ e^''v\r,\''e-^Pdfigdt 

>2 [ e2^*(7-(5)||?7||^i dt - 2sup||7;||i2 [ e'^^^Ml^ dt 
Jo ^ t Jo 

2 / e2^* 

> / e^^\2j-26 - 4C,2sup||v||i2)||r/||^i dt 

Jo t A 


as there exists Cg > such that \\r]\\n < Cs\\'r]\\f^i by the Sobolev inequality. (Covariant Sobolev 
spaces and extensions of standard inequalities in these were discussed in the appendix in j3]). By 
the energy non-increase, 

sup \\v\\l2 < V2(V - ttN) < \/2(Vo - vrA^) = (3.27) 


and so for y/2eo < "^lyj^^ (which is implied by the assumption on eg), the final term on the right 
hand side of (|3.26p is bounded below as 


> 6 / e2^*||r/||^i dt. 
Jo ^ 

Hence from (j3.23p 

sup e'''Mt)\\l. + S f e''%\\l. dt< Uml^- (3.28) 

tGfO.Tl JO ^ 

Using this and (|3.27p (together with Holder's inequality and the same Sobolev imbedding as above) 
we can bound the final term on the right hand side of (j3.24p : 

8 /7 e'''v\rj\\-'^d^,gdt<8 r [ e'^'M,. C^Ml. dt < ^^^\m\\h < imWh 
Jo Jy, Jo Jt. Ad 

for ^/2e^ < Using this and in we have 



5 e^''\\v\\l,dt <2 e^''{Q^{v)-6\\v\\l,)dt < [WmWh + \Hmh) (3-29) 

which proves the corollary. □ 
In terms oi v = B — ^(1 — and r/ = Ba^ = ^{Di + iD2)^, the equations (11.60 read (in 

gauge invariant form): 

dAi dAo . dv 

dA2_dAo^ (3.30) 

We will now show that the equations p.30p together with the estimates p.25p and also (jl.lOp 
for Aq imply the convergence as t — )• oo in gf 1 1— ). (A, $)(t) at an exponential rate. 

Lemma 3.4 In terms ofr],v above, 

\\A\\L2 + \mL2<cm\HX + \Hm) (3-31) 

which together with the estimate (j3.25p completes the proof of the theorem in the case |S| > 4ttN . 

Proof Unless specified otherwise the generic constants c, c' depend only on S. 

;i _ |$|2)2 = _ 2||$||^2 + ll^ll^^ > |s| - 2||^>||^4|S|5 + ll^ll^^ > i||$||^4 - 

(as ||$||22 < ||$||2^|S|3 < i||$||^4 + |S|). Therefore, 

V(A,$) > c(||L»A^|li2 + ll^lli*) -c' > c||$||^^-c' 

and using the energy non- increase (shown at the end of the proof of the existence theorem 12. ip this 

||$||i4 <c(Vo,S) (3.32) 

where Vq is the initial energy. 

The last equation in (j3.30p implies 

\\Donh < cdir/ii^^ + ib^iii^) 
< c{M\i,^ + \mUv\\%) 

<c(Vo,S)(||r7||^^ + ||H|^,) 

Therefore by (13:25]) . 

Do^Wlie'^^^ dt <oo (3.33) 

which implies the exponential decay of t i— )• Dq^. However, as seen in the previous section, the 
equation for <I> implies in the Coulomb gauge 

and then the Calderon-Zygmund inequahty, given the condition f Aodfig = 0, imphes 

\\Ao\\H2,r, < c{^,p)\\{i<^>,Do<^)\\LP. 

For p < 2, 

\\{i<P,Do<^)\\Lv <\m (3.34) 

<c(Vo,S)||L>o$||i2. (3.35) 

Therefore, together with (|1.10p 

Po|ll/2,p + lld^ollia < cpo^llia. (3.36) 



/ (Pollia.p + \\dAo\\l2)e^^' dt < oo (3.37) 

by ()3.33p . By the equation for $ in ()3.30p we also have 

II— $11^2 < po^>||i2+c(||??||^^ + ||?;||j^i) (3.38) 

< c(Vo, S) (Po^IIl^ + ll^ll^i + MhA (3.39) 

and so 

Similarly for A, we obtain 

I ^ ^||2 „2St 



dt < oo . (3.41) 

These estimates immediately imply ()3.3ip by the equations (j3.30p and then by ()3.25p we conclude 

^(A, <!>) = / on [0,oo)xS 

where / G and ||/||j;^2 (t)e^'^*(it < oo, and thus 

(Aoo, $oo) = limt^^{A{t), ^>(t)) 

exists in with exponential rate of convergence. 

From these estimates and the above we have that ^> G and so is in for all p < oo and 
as i — ^ oo, <I> <^oo strongly in L^, Da^ -Daoo^oo weakly in L^, hence ^Da^ — > ^I'oo-Daoo^oo 
weakly in L^; in addition, as <I> is bounded in L"? for any q < oo (from <I> G H^), the product 
^Da^ is bounded in L'' for every q < 2. This implies weak convergence in L'' for g < 2 so that 
div{i^, Da^) converges weakly in H~^''^; by Calderon-Zygmund then Aq converges to a limit Aqoo 
weakly in H^'*^ for each q < 2 and so strongly in every for p < oo. Clearly by (13.370 this limit 
is zero. 

This completes the proof of the theorem in the case |S| > 47rA^. □ 

3.2 Proof in the case |S| < iirN 

In this case a different approach is needed because the condition ||^||^2 ^ ^ > which held and 
was used in the previous section is no longer valid (as is obvious from the explicit form of the 
lower bound in ()3.22p ). Indeed, as first observed in [2j, in this case there are no solutions to the 
Bogomolny equations ()1.5p . To see this fact, integrate the equation v = over S to deduce 

2nN - ^ + = (3.42) 

which is an impossibility when 47rA^| > We will show instead that the minimizers of V have 
u = 1^ = constant and H'l'lli^a = 0, i.e. the field <I> is identically zero and Da = V — iA is a 
constant curvature connection on L. Notice that the value of the constant / is fixed as 

/ = 27riV - > 0, (3.43) 
since vd^g = / + |||'&||^2 by integration of (I3.17p . Now define d> ^hy d? = \T,\/ v'^dfig, then 
2^y w^d/Xg^ ^ |$]| 2 = d J w^d/Xp + -Li = 2 J vdfig + J [Vdv — —=)'^d^g 

= 2{i + ^\ml,) + \\{Vdv-^)\\l,. 

Referring to (jl.4p . this implies that the energy V can be written 

V(A, $) = 2\\dAnh + ^ [2(/ + ^Wnh) + U^v - + ^N. (3.44) 

From this it follows, since / > 0, that V > ttN + ^j^y that this lower bound is achieved 
with $ identically zero and B = = constant, so that v = = ^ = constant. Thus although 
the Bogomolny bound vrA^ is not itself achieved in the case 4:7rN > |S| we have identified the 
greatest lower bound and shown that it is achieved with the "pure magnetic" constant curvature 

Lemma 3.5 For \T,\ < inN 

■ .r (47rA^-|S|)2 

mm V = vriV + , ' 

m{L®n^) 8\T,\ 

and this minimum value is achieved by taking the Higgs field ^ identically zero and D to be a 
constant curvature connection. 

To analyze the gradient flow in this case it is useful to introduce the variable y = v — I (in place 
of v), so that (|3.18p are replaced by 

Doll - AdAie-^PdAV) + {l + |$P)r? = -yrj, 

(^_A+|cI>|2)y = -|<I>|2/_4|^f 

Do<S> - Ae^'^PdAdA'^ + l<^ = 

|--A + |cI>|2)y = -|<I>|2/_4|,?|2, (3.45) 

Corresponding to (j3.23p we have the fohowing integral inequahties for solutions of (j3.45p 

y(r)||i,+2^ e''''[Q^iy)-5\\y\\l,)dt < ||y(0)||i. - 2 ^ ^ e^^d^l'/ + 4|r?|2)e-2^d/i,dt, 

2^^||<I>(r)||22 +2 j"^ e''^'(\\dAnl2 + {l-6)\ml2)dt < \m0)\\l,-2 j"" j e^^' y\^^ e-^P d^igdt. 


Notice that since I > 0, there is a natural mechanism forcing $ to converge to zero at an 
exponential rate. However this fact necessitates modification of the arguments based on the lower 
bounds in lemma [321 which are dependent upon ||$||2,2 > m > 0. The presence of / > in the 
integrals in the first and third identities means that for 5 < | there exists 7 = 7(Vo, /) such that 

QAMv) + {l-m\\h>7h\\li^, and (3.47) 

WdAnh + (/ - S)\ml, > mli^- (3-48) 

For the middle identity, recall that ify = -^Jj.y then Poincare's inequality says that 

/ e-^P Yl \9ky\^dfig > Cp\\y - y\\l, > Cp{\\y\\l, - \\y\\l,), 


SO that there exists 7' > such that 

But also I + y = 1^/2^ = ^+ 2|b| H^Hi^' ^° ^^sX y = 2l^ll'^lli2- Therefore, overall we have the 
following set of controlling inequalities: 

e^'^MT)\\l,+2 r e^^'^Ml^dt < ||r?(0)||i, - 2 T / e^M^pe-^Pd/i.dt, 


' ^ Jo JT. 

2^^"y(r)||i.+2^%2^*((7'-<5)||y||^.-yily|li.)dt < Il2/(0)||? 

-2 r / e%(|$p/ + 4|7?|2)e-2^d/i3dt, 

Jo JT. 

rp rp 

^T)\\l, + 2 [ e^^'^m^dt <\M0)\\l,-2[ [ e^^'y\<^\^e-^fdfigdt 
Jo ^ Jo Jt, 

Add the three inequalities, and use the fact that ||y||^2 = 1^1 |yP = ^y||^||^2 to absorb the negative 
term in the integral in the second inequality by the integral in the third one. Next bound the 
nonlinear terms on the right hand sides in the same way as in (j3.26p - (|3.25|) . and conclude that as 
long as 112/(0)11^2 is initially small then 



e^'^MT)\\l, < \m\\h, 

e''''UT)\\h < wvmh, 
e^'^mmi, < wmwh, 

holds for all T > 0, and 




y\\m + \\nli^ + \\v\\li^ 

dt < oo 


from which convergence can be deduced as in the previous section: (j3.33p and (j3.37p hold as a 
consequence of (j3.49p . and hence by the first two equations of (j3.30p the estimate (j3.4ip also holds, 
and so A(t) converges to a limit A^o at an exponential rate. The proof of the main theorem is now 


[1] E. Bogomolny, Stability of Classical Solutions, Soviet Journal of Nulclear Physics 24 861-870 (1976). 

[2] S. Bradlow, Vortices in holomorphic line bundles and closed Kaehler manifolds, Commun. Math. Phys. 
118 1-17 (1990). 

[3] S. Demoulini and D. Stuart, Adiabatic Limit and the Slow Motion of Vortices in a Chern-Simons- 
Schroedinger System, Commun. Math. Phys. 290, 597-632 (2009). 

[4] S. Demoulini and D. Stuart, Gradient flow of the superconducting Ginzburg-Landau functional on the 
plane, Commun. Anal. Ceom. 5(1) 121 - 198 (1997). 

[5] S. Demoulini, Global existence for a nonlinear Schrodinger-Chern-Simons system on a surface, Ann. 
Inst. H. Poincare Anal. Non Lineaire 24(2) 207-225 (2007). 

[6] E. Feireisl and P. Takac, Long-time stabilization of solutions to the Ginzburg-Landau equations of 
superconductivity Monatsh. Math. 133, 3 197-221 (2001). 

[7] A. Hassell The Yang-Mills-Higgs heat flow on J. Funct. Anal. Ill (2) 431-448 (1993). 

[8] A. JafFe and C. Taubes, Vortices and Monopoles, Birkhauser, Boston, 1982. 

[9] J. Jost, Riemannian geometry and geometric analysis, Springer- Verlag 1988. 

[10] M. Noguchi, Yang-Mills-Higgs theory on a compact Riemann surface J. Math. Phys. 28, 2343 (1987). 

[11] L. Simon Asymptotics for a class of nonlinear evolution equations, with applications to geometric 
problems Ann. of Math. 118 (3) 525-571 (1983). 

[12] M. Taylor, Partial Differential Equations I-III, Applied Mathematical Sciences, vol 117, Springer- Verlag